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Abstract: The low multipole anomalies of the Cosmic Microwave Background has re- 
ceived much attention during the last few years. It is still not ascertained whether these 
anomalies are indeed primordial or the result of systematics or foregrounds. An example 
of a foreground, which could generate some non-Gaussian and statistically anisotropic fea- 
tures at low multipole range, is the very symmetric Kuiper Belt in the outer solar system. 
In this paper, expanding upon the methods presented in [Q], we investigate the contribu- 
tions from the Kuiper Belt objects (KBO) to the WMAP ILC 7 map, whereby we can 
minimize the contrast in power between even and odd multipoles in the CMB, discussed 
mil- 

We submit our KBO de-correlated CMB signal to several tests, to analyze its validity, and 
find that incorporation of the KBO emission can resolve the quadrupole-octupole alignment 
and parity asymmetry problems, provided that the KBO signals has a non-cosmological 
dipole modulation, associated with the statistical anisotropy of the ILC 7 map. Addi- 
tionally, we show that the amplitude of the dipole modulation, within a 2a interval, is in 
agreement with the corresponding amplitudes, discussed in j|. 
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1. Introduction 



The departure of the CMB from statistical isotropy and homogeneity has attracted very 
serious attention since the first publications of the WMAP data || ^, \j§ . Investigations of 
the low multipole anomalies of the WMAP co-added map and the ILC whole sky map can 
clarify possible sources of "contamination" g g |, [0], 0, g p|, 0, [D| 0, ||, |||, |0 



21]. The low multipole non-Gaussian features and departure from statistical anisotropy of 
the CMB can be related to the foreground residuals |l|, [2^, 23, |4|], systematic effects p5| . 
or it could even have a primordial origin, for instance by a violation of the Copernican 
principle [26|, a primordial magnetic field [27, pH , or by a non-trivial topology of the 
Universe |2^, |30|, ^TJ . Usually, discussing uncounted residuals of the foregrounds as a source 
of low multipole anomalies, one can assume that most of those residuals are associated with 



the Galactic plane, while for instance, the observed quadrupole -octupole alignment [16] is 
connected with the ecliptic plane. Spergel et al. |3^] have pointed out, that low multipole 
anomalies, like a power asymmetry and statistical anisotropy of the CMB signal, correlates 
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with the ecliptic North and South poles, and more generally, reflect the morphology of the 
map of number of observations (MNO). Indeed, [33, 34] building on the work of 35] and 



1 36] have identified a quadrupolar power asymmetry, confirmed by the WMAP team 
Most likely, this anomaly indicates an influence of the beam asymmetry on the CMB signal, 
discussed in pq] . However, as it was pointed out by |HJ], in spite of significant distortions 
of the phases of the CMB signal, an anisotropic beam could only change the beam window 
function with less then 0.6% at I 3> 400, and it is insignificant for estimation of the CMB 
power spectrum at the low multipole range I <C 200, with changes at the level of < 0.1%. 
Most of these tests are based on the CMB maps, meaning that a high quality of data 
reduction is essential. In contrast, as recently proposed in J2|, ||, the parity test requires only 
the CMB power spectrum, which is usually estimated with significantly better accuracy 
than the CMB map |gcj. 

The idea of the parity test is based on the analysis of the ratio g(l) of the powers, stored 
in even and odd multipoles. For a statistically homogeneous and isotropic random CMB 
signal, <?(/) should have no preferences in respect to I = even or I = odd (see [41, 42] for 
discussions). 

In this paper we would like to focus on an extension of the parity test, quadrupole-octupole 
alignment and lack of angular correlations at 9 > 60° 1 for the whole sky CMB map, taking 
under consideration the possible contributions from "new foregrounds" like the Kuiper Belt 
objects KBO, (see for review [43|] ), counting them as a source of contamination of the CMB 
Q. Under the standard assumption that the KBOs is localized mainly in the Ecliptic plane, 
this foreground is potentially "dangerous" for the quadrupole-octupole alignment. More 
importantly, due to the high degree of symmetry, the morphology of the KBO signal could 
be remarkably similar to the morphology of the MNO, for instance at the V-band, and at 
least at the low multipole range I < 20 — 30, where the power is comparable to the CMB 
power spectrum. 

The critical point of our analysis is the assumption, that the dipole modulation of the CMB 
Q, is associated with an unknown systematic effect. Thus it affects not only the measured 
primordial CMB signal, but all the measurements of the foregrounds, including the KBO, 
as well. 

We estimate the temperature change which the model-KBOs contribute to the CMB. We 
show that unlike the model of Q, de-correlation of the ILC 7 and KBO leads to a decrease 
of the power of even and odd harmonics, improving the shape of the parity parameter g(l). 
Moreover, we will show that a KBO-ILC 7 de-correlation could change the significance of 
the quadrupole-octupole alignment, down to the level of spontaneous (chance) correlations. 
Thus, if the KBOs are responsible for the parity asymmetry of the CMB, one should be 
cautious of the alignment of quadrupole and octupole components of the CMB. 
The outline of the paper is the following. In section 2 we introduce the basic characteristics 
of the parity asymmetry. In section 3, we discuss the model of the Kuiper Belt. Then, in 
section 4, we present a method of cross-correlations of the KBO foreground and the WMAP 
ILC 7 map at low multipole range. Section 5 is devoted to the analysis of general properties 

1 As it was shown in the anomaly of the correlation function can be explained as manifestation of 
the parity asymmetry of the CMB power 
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of the cross-correlations between ILC 7 and KBO signal. In section 6 and 7 we investigates 
the distortion of even and odd component of the ILC 7 map. Section 8 is devoted to the 
cleaning of the ILC 7 map by de-correlation with KBO emission. Additionally we look 
at the implications for the parity asymmetry and the quadrupole-octupole alignment, and 
test various properties of the KBO filtered signal in Section 9. Finally, in section 10, we 
summarize our results and draw our conclusions. 



2. Odd multipole preference of the power spectrum 

The temperature fluctuations of the CMB can be decomposed into spherical harmonics in 
the standard way: 

oo I 

AT(0, <^>) = E E a hrnYi, m {e, 4>), (2.1) 

Z=l m=—l 

where a^ m is the coefficient of decomposition: a^ m = \ai iJn \ exp(«<^ )m ), with <j>i >m as the 
phase. 

Under the assumption of total Gaussian randomness, as predicted by the simplest infla- 
tionary model, the amplitudes \ai >m \ are distributed according to Rayleigh's Probability 
Density Function (PDF) and the phases of a^ m are supposed to be evenly distributed in 
the range [0, 2vr] @. 

For any signal T(n) defined on the sphere, one can extract a symmetric T + (n) = T + (— n) 
and an anti-symmetric T~(n) = — T~(— n) component, where 

oo I 

T± ( A ) = E E ^, m r ± (/)y,, m (n), (2.2) 

1=1 m=—l 

and r+(Z) = cos 2 (f ), T~(l) = sin 2 (f ), Y l>m (A) = Y l>m (-h) . 

For the concordance ACDM cosmological model with initial Gaussian adiabatic perturba- 
tions, we do not expect any features distinct between even and odd multipoles. However, 
there have been reported power contrast between even and odd multipoles of WMAP TT 
power spectrum ||, 37, 45]. The corresponding estimator for "even and odd" asym- 
metry of the CMB power spectrum is given by [§]: 

V± n(n + l)C(n)T + (n) 

g(l) = V — — (2.3) 

T!n min n{n+l)C{n)T-(n) 

where C(n) is the usual power spectrum. At lowest multipoles (2 < / < 23), there is 
odd multipole preference (i.e. a power excess in odd and a deficit in even multipoles) as 
discussed in 10, 0, Eol M. 



3. KBO's as a 'new foreground' 

The contribution from the residuals of the foregrounds into parity asymmetry of the CMB 
was widely discussed in [l], ^, E3]. Unlike the first two papers, the last one claims that 
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Figure 1: Left column. Model for the KBO emissivity distribution (arbitrary units) in Ecliptic 
coordinates for H = 15° (top), for H = 30° (the middle panel) and H = 70° (the bottom panel). 
Right column show the KBO emissivity in the Galactic coordinates. All maps carry equal power. 



uncounted foregrounds, like emissivity of dust from the Kuiper Belt, could explain the 
detected parity asymmetry of the CMB. The main idea of the method proposed in Q 
is that a very symmetric foreground, when taken into account and subtracted from the 
WMAP whole-sky ILC map, could amplify the CMB power stored in even multipoles, 
increasing the parity parameter g(l) (i.e. mitigating the observed parity asymmetry) for 
the multipole range 2 < I < 23. An illustration of the idea is shown in Fig. [l], where the 
models of the KBO emissivity are shown in Galactic and Ecliptic coordinates for different 
angular heights of the KBO H = 15°, 30° and 70°. 

The blackbody-like radiation from the KBO affects the intensity of the microwave sky as 
follows: 



I{y, f ) = {B{y, Tkbo) - B{y, T C mb)) r(f ) 



(3.1) 



where v is the observation frequency and B{u,T) is the blackbody radiation spectrum at 
temperature T, and r(r) is the optical depth of KBOs at the sky direction r. The second 



term on the right hand side of Eq.(3.1) arise from the occultation of CMB photons by the 
KBO. The FIRAS data imposed the most stringent constraints on the sky-averaged optical 
depth of KBOs, which is r <, 3 x 1(T 7 Q g|, @ HI- 

The heating is mainly due to the radiation from the Sun. Therefore, it is possible to 
calculate the temperature of these objects at the equilibrium, assuming that it arises from 
the conversion of the solar radiation absorbed by the object into microwave emission. At 
a distance of 40AU, where KBOs are most densely populated, we find an equilibrium 
temperatures of ~ 43. 7K 50|. In Table p], we show the frequency spectrum of KBO 
divided by the CMB anisotropy spectrum f{v). 
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As evident, the frequency spectrum 
of KBOs does not vary significantly even 
between the both ends of the WMAP 
observation frequency: 

(fM - fM)/fM & 0.22,(3.2) 



band 


K 


Ka 


Q 


V 


W 


m 


41.48 


42.07 


42.67 


44.8 


50.42 



Table 1: Frequency spectrum of KBOs at WMAP 
frequencies (in GHz) normalized to CMB anisotropy 
spectrum. 



Therefore, the KBO emission may be confused with intrinsic CMB anisotropy, even when 
the contribution to the microwave sky emission is sizable. For the allowed values of the 
optical depth, r <, 3 x 10~ 7 , the KBOs may have an averaged effect on the CMB data as 
big as ~ 15 /xK. Provided the KBOs have certain large-scale patterns, KBOs may have an 
effect large enough to be the culprits of reported large-scale CMB anomalies [Q, ^ [14], 16 

i7|, 0, ig, g§ gg, n, g, n n g, [55]. 



3.1 Angular distribution of the KBO emissivity 

To investigate the emissivity of Kuiper Belt, we adopt the models by yj, where the Kuiper 
Belt covers a symmetrical band in the Ecliptic plane of constant disk height H, with 
uniform temperature distribution. Following [|J, we pick three models, defined by the angle 
±Hkbo/2 from the Ecliptic plane, where we choose Hkbo = 15°, 30° and 70° respectively. 
For the three values of Hkbo, we set the temperature inside the Kuiper Belt, so that it 
agrees with the upper bound of 15/uK for the entire sky, as calculated previously in the 
article. In practice, we do this by calculating the respective fractions of the total sky area, 
for each value of Hkbo, an d finding the required temperature inside the band from this 
value: 



- kbo\ 



Be 



e 



(3.3) 



where B is a normalization constant and O(x) is the Heaviside- function. In the discussion 
which follows, we will use the average temperature of the KBO emission as normalization 
of B : 



Tkbo = h 



2tt 



d9 sin 9Th 



kbo( 



„ (H 
B sin — 
V2 



(3.4) 



Now, we decompose the KBO signal into spherical harmonics and get the following coeffi- 
cients of decomposition: 



fl 



Lm 



21 + 1 

47T 



BT+(l)S m , 



(ir+H)/2 



(*-ff)/2 



d9 sin 6* P/ (cos ( 



(3.5) 



For H 7r/2 in Eq.( |3.5|) we can use the following Taylor series representation: Pi(x) ~ 
P/(0) + ±P{'(0)x 2 , where the two primes denotes the second derivative. The first derivative 
vanishes for even I = 2n, due to the symmetry of the model ( dP ^ U=o = 0). This gives 



fi(0) 



r(^i)r(^: 



-1 



,P/'(0) = -/(/ + l)PK0) 



(3.6) 
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Here T(x) is the Gamma-function. For 1(1 + 1)H 2 /24 <C 1, from Eq^SJj-SJ)) one can get: 



/,.„. > M /^i? S in(|)P,(0)r+(/) 



1 



^ + 1) ^2 

24 



(3.7) 



Thus, for /« 35 (^) 
representation: 



flfi 



the coefficients of the expansion fa q have the following analytical 
2/ + I.„ . „(2n-l)!! 



7T 



-T KBO (-l) 



2"n! 



£ = 2n 



(3.8) 



In order to test the influence of our models on the lowest multipoles, we now find the power 
spectrum for each value of Hkbo an d compare them with the ILC power spectrum. The 
results are presented in Fig. |2[ It is clear, that for all 3 values of Hkbo-, the power of the 




Figure 2: Power spectrum D(l) — K^ 1 ) C(l) (in units of mK) for the three KBO models, with 
Hkbo = 15° (black) 30° (red) and 70° (blue), compared to the ILC power spectrum (black line). 

KBO quadrupole is strong enough to affect the ILC quadrupole. For higher multipoles, 
Hrbo = 15° can affect the even multipoles up to I = 19. 

4. CMB-KBO cross-correlation 

In this section we will discuss general properties of the KBO foreground, basing our anal- 
ysis on the symmetry in the Ecliptic system of coordinates. Firstly, we will assume that 
the emissivity of the foreground S(9 e ,4> e ) is very symmetric with respect to the plane 
6 e = ir/2 in Ecliptic coordinates, and do not depend on the azimuthal coordinate (j) e : 
S(6 e ,4> e ) = S(6 e ). This assumption, leads to the following properties of the foreground: 
fi,m = /(0^m,o^/,2n) n = 1, 2..., clearly demonstrating that only even multipoles can be 
affected by the KBO-foreground. 

At the same time, at the level of 3a, the ILC 7 map has dipole modulation, which will 
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provide a coupling between the even and odd components of the KBO emission, depending 
on the orientation of the dipole term and the corresponding amplitude of modulation. Here 
we assume that this kind of modulation has a non-cosmological (systematic) origin, leading 
to a distortion of the CMB and the foregrounds as well. 

To assess the problem of a possible contamination of the ILC by the KBO-foreground, we 
will use the model of the CMB signal for a given direction on the sky n: 

T c (n) = T ilc (n) + F { Xkbo (n) + e(n)} 

= T ilc (n) + [1 + T(qn)] X kbo(n) (4.1) 

where T c (n) and T^ c (n) are the temperature of the intrinsic CMB signal and the ILC, 

Xfcbo(n) oc Tfcfe (n) corresponds to the residuals of the KBO emissivity, F{..} denote the 

functional of the variables Xkbo( n ) an d s(n), T(qn) is the function of dipole modulation, 

q is the direction of the dipole modulation, and e(n) = Y(qn)x£;6o( n )- 

In the forthcoming discussion we will assume that T(qn) = A(qn), where A = const and 

qn is a dot product of two vectors q and n, where q points in the direction of the dipole 

modulation of the CMB. Thus, in our model, the second term in Eq. fl4.lD , the residuals of 

the KBO emissivity, is symmetric with respect to inversion: Xkbo( n ) = Xkbo{~ n )i while the 

last term, dependent on the dipole, is anti-symmetric: e(n) = — e(— n). 

In the multipole domain these two component contribute to the coefficients of the spherical 

harmonic decomposition additively: 

Q,m = ai,m ~ f(0Xi,m r+ (0 ~ P(0 £ J,m r ~(0 = a l,m ~ ^l,m (4.2) 



where ci >m is the primordial coefficient of decomposition, are given by Eq. (|2.2| ), q(l) and 
p(l) are weighting coefficients, and e/ m can be expressed as 



£ lm= X/ ( _1 ) m ^lm'/i"m") 

(4.3) 



m'=-ll"=0m"=-l" 




3(21" + 1) (21 + 1) I" I I 

4^ V / 

Here b\ m relates to the parameters A and n by the relation b\ m = ^-l^(n). 

As stated earlier, for the given model of the KBO we have fi :Tn = f(l)T + (l)5 m fi, and thus, 

in Ecliptic coordinates, Eq.(fO|) can be simplified as follows: 



e lm = > ' > (-ir^(l")b lml f v r' 3(2? " + 1)(2Z + 1 




47T 

(4.4) 



5. General remarks about KBO -intrinsic CMB cross-correlation 

In this section we would like to show, that the model can work if, and only if, the primordial 
CMB signal just by chance has a very strong cross-correlations with uncounted foregrounds. 
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Let us discuss this issue in more detail. Following Eq.( [4.2[) , we assume that the ILC 7 signal 
can be represented in the following way 



ai,m = ci >m + ILi )m (5.1) 

Let us define the coefficients of cross-correlations between the ILC map and the foregrounds, 
and the primordial CMB and the foregrounds as follows: 

a l,m^*i m + a l tT rfll,m 

K{1) 



i > 

2\2 



2(E m l«i,m| 2 Em' \ U l,m'\ 2 ) 
Em c l,m^i m + c lrrJ^l,m 

k(1) = l - j - (5.2) 



2(£Jci,J 2 £ m ,|n iim ,| 2 )* 

Combining Eq.fl5.1|) and Eq.( |5.2[ ) we can get the following equation for the <r(Z) -parameter: 



12x1 /V- U. 12X1 



t(0 = Jrro (&fep)'-* ( ' , (fefel (5 - 3) 



So, as one can see from Eq.(|5.3|), the coefficient play a role of effective coefficient of 
cross-correlations between ILC, intrinsic CMB and the foreground residuals. From Eq.(|5.1|) 
and Eq.( |5.3[ ) one can define the power spectrum of the primordial CMB as: 

C P (l) = E = V{l)C ilc (l) (5.4) 

m 

where Cu c (l) is the power of the ILC 7 map, and 

is the factor of modulation of the ILC 7 power. Now, by using Eq.(|5.4|) we can estimate 
the parity parameter for the primordial CMB as 



m = Y!n=2<n + l)V{n)T + {n)C llc {n) 
"" ! ~ Yln=2n(n + l)V(n)T-{n)C Uc {n) 



In order to increase the contribution of even multipoles in parity parameter g p (l), the func- 
tions V+(l) = V(l)T + (l) and V~{1) = V(l)r~(l) should satisfy the following conditions: 
V + (l) 3> V~(l). That means that there are only two variants: 1 or \K~\ — > 1, 

where the + and — in the superscript symbolizes that the n or the K is drawn from either 
V + or V~ respectively. The first case (with — > 1 and symmetric foreground) was 
discussed in pj, and it requires a significant coupling between the primordial CMB and 
KBO-foreground at the level \k + \ ~ 0.7 — 1. 

Due to the very clear established non-Gaussian properties of the KBO-foreground, a high 
correlation between it and the primordial CMB temperature anisotropy, is in contradiction 
with the assumption about a statistical isotropic and Gaussian intrinsic CMB signal. Thus, 
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the improvement of the parity parameter by an increase of \k + \ (hereafter Model 1 ) will 
have lost its theoretical basis. 

In contrast, for |k ± | <C 1, the only way to increase g p (l) from Eq.( |5.6| ), keeping the as- 
sumption about Gaussianity and statistical isotropy of the intrinsic CMB signal, is to get 
\K~\ —7- 1 and \K + \ < 1 (Model 2). That is, a high correlation between the ILC signal and 
the KBO-foreground for odd multipoles, and at the same time, a relatively low correlation 
between the ILC signal and the KBO-foreground for even multipoles. 

6. Distortion of even multipoles. Estimation of K + in Model 1 

In this section we will estimate the coefficient of cross-correlation K + for the ILC 7 and 
the KBO foreground for even multipoles, taking into account Eq.( |5.1| - |5T2"| ), Based on the 
symmetry of the KBO, we will perform our analysis in the Ecliptic system of coordinates. 
From Eq.( |5.1| ) we get: 

ci,o = %o- s(l)fl, ci,m = ai,m, ,m>l, 
C p (l) = C ilc {l) + (? 2 /f - Kawfi) 

= C Uc (l) + C f (l){,\l)-2,(l) a -^j (6.1) 

where C/(Z) is the power spectrum for the foreground data. We can now express the last 
term in Eq.(^) as: 

= C p (l) - Cg c (l) = Cf(l) ^ ^ 

Cilc(l) Ca c (l) 

where 

w(Z) = ? 2 (Z)-2 ? (Z)^. (6.3) 
Jl 

Let us define the coefficient of cross-correlations for ILC and the foreground as 

0(0 = M™M.oW*'.o) ( 6 .4) 

(2Z + 1)IC| C (0 

where ipi o is the phase of the m = modes for the ILC and = ^(l + cos(nl)) is defined 
as phase for the foreground. In Fig. ||| we plot the coefficient Q(l) versus I in Ecliptic 



coordinates from Eq.(6.4). As one can see from this figure, the range 25 < I < 40 reveals 
only negative correlations, which looks quite strange for a random Gaussian process. This 
range of multipoles requires special investigation, which we will present in a separate paper. 

7. Distortion of odd multipoles. Estimation of K~ in Model 2 

The contribution from the KBO emissivity into the coefficient of cross-correlations K~ = 
K(l)T-(l) is given by Eq.Q and Eq.|T 



Em a l,m £ *m + a Lm £ l:> 

K-(l) = l - ^, (7.1) 

2 (Em Km| 2 Em' \ £ l,m'\ 2 )'' 
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Figure 3: The coefficient Q(l) versus I in Ecliptic coordinates. 
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Table 2: as jTn coefficients of the ILC 7 octupole in the Ecliptic coordinates. 



where 



e i>m = AcosQ[ai(l,m)fi +1 , m + Pi(l,m)fi- ljm ] 
+ (A/2) sine e^[a 2 (/,m)/^ 1)m+1 - 
- /5 2 (/,m)/ J+ i )m+ i] 
+ (A/2) sin 6 e-^ia^m)/^^ - 

~ ^3(/,"l)/i-l,m-l], 

(7.2) 
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and 



ai(l, m) 
a 2 (Z,m) 
a 3 (l,m) 



(l + m + l)(l-m+l) 

(2Z + l)(2/ + 3) 
(/ + m)(l — m) \ 1 I 2 



1/2 



(21 + l)(2l - 1) 
(I — m)(l — m — 1) \ 



(21 + l)(2l - 1) 
(/ + m + l)(/ + m + 2)^ 1/2 



(2/ + l)(2/ + 3) 
(I - m + - m + 2)\ 1/2 
(2Z + 1)(2Z + 3) 



^ m)= V (2/ + l)(2/-l) J • (7 - 3) 

Here and $ are the coordinates of the dipole modulation in the Ecliptic system of 
coordinates, and fi j7n = fiT + (l)5 m fi is the KBO foreground. Due to the high symmetry of 
the KBO foreground, only the even multipoles I = 2n, n = 1, 2, .. are presented in Eq.(7.2- 



7.3| ). That means that £; im -term in Eq. (|7.2| ) has only three non- vanishing component with 



odd I = 2n + 1 and m = 0, ±1. Namely, 



£2th-i,o = Acos@[ai(2n + l,0)/ 2n+2 

+ A(2n + l,0)/ 2n ], 
£2^i,i = (A/2) sin 9 e^[a 3 (2n + l,l)f 2n+2 

- /? 3 (2n + l,l)/ 2n ], 
£ 2n+1 ,_i = {A/2) sin 9 e^[a 2 (2n + 1, -l)/ 2n 

- /? 2 (2n + l,-l)/ 2n+2 ] (7.4) 



where: Z = 2n + 1, and 



.2,,-^ „2„ ^ ^ + 1) 



a£G,l)=$(Z,-l) 



(2Z + l)(2/ + 3)' 
^,l) = ^,-l)= (2Z + 1)(2/ _ ir 

= (2/ + lK2! 2 + 3) '^^ 0) = (2/ + 1K2/-1) - ™ 

7.1 Particular model of dipole modulation. Maximization of the octupole com- 
ponent of cross-correlation 



As it is follows from Eq.( [7.l|) and Eq.(7.2), the coefficient of cross-correlations depends on 
orientation of the dipole in the Ecliptic coordinates (0,<3?). Let us discuss one particular 
choice of and <£, which maximize the coefficient of cross-correlations K~(l = 3). As 
we have pointed out in introduction, our main goal is to show that KBO foreground can 
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resolve the parity problem at the range of multipoles 2 < I < 23. However, this model 
can potentially explain the quadrupole-octupole alignment also, due to coupling between 
even and odd components of £i m , given by Eq.([7.4|). In particular, the octupole compo- 



nent of £/ jm , according to Eq^T^) depends on the linear combination of the quadrupole 
and I = 4-components of // jm . Due to the azimuthal symmetry of the KBO emissivity, 
these components have opposite phases. Thus, the coefficient of the KBO quadrupole and 
octupole cross-correlation depends only on O and <F Since KBO quadrupole and octupole 
components could correlate with the corresponding components of the ILC, maximizing 
correlations between the ILC octupole and £3, m , will provide the optimal way for a dis- 
alignment of the intrinsic quadrupole and octupole. 

In Table § we show all the component of the ILC octupole in ecliptic coordinates. As one 
can see from this Table, the phase of (3, 1) component of the octupole is 923,1 = 1.227, 
which means that in order to maximize the K~{1 = 3)-coefficient, the azimuthal angle $ 
should satisfy the following equation: = — <ps ±. Then, the coefficient of cross-correlations 
K ~ (/ = 3) is given by: 

CI3 ob cos O + 2|a3 i|csin O 



7C(l = 3) [b 2 cos 2 + 2c 2 sin 2 9] 



a 



2 iJ. 1 ,l„2|„_ . |2\ 2 



3.0 



V +4c^|a 3 ,ir \ 2 cos(0-r?) 



(7.6) 



\ 7C(/ = 3) J cos 2 + 2c 2 sin 2 6 

where we have used the following definitions: £3^ = b cos© and (£3,1 1 = csinO (see 
Eq.(^J)), C(l = 3) is the power of the ILC octupole and 



n = tan- ( M 1 . (7.7) 

V «3,0» / 

The nominator in Eq.( |7.6|) has a point of maxima at O = 77, and it vanishes at O = rj + n/2. 
As it follows from Table ^, the amplitudes of the octupole in the Ecliptic coordinates are at 
the same order of magnitude, while the parameters b and c depend on the power spectrum 
of the KBO angular anisotropy. This is why we will focus on two asymptotics, |6| S> |c|, and 
\b\ <C |c|, in order to investigate the dependency of the cross-correlations on the particular 
choice of these parameters. In the case \b\ S> |c|, all the power of the £; jrrt -signal would be 
concentrated at m = mode, and 

, Q3 '° ^-0.65, ©«-. 

sTm = 3) 2' 

In the opposite case, when |6| <C |c|, the coefficient of cross-correlation is given by 

In Fig. ^ we show the dependency of i^ z ~ 3 (0) for different KBO models with H = 
15°, 30°, 70°. As is seen, for the KBO models with H = 15° and H = 30°, the point 
of maxima of Kr s (0) is close to O ~ tt, while for H = 70° it is close to 7770 — 138° from 
Eq.([7.7|). Next Fig. ^ illustrate the cross-correlation coefficient for KBO H = 15°, H = 30° 
with O = 7r. 



^/= 3 (Q) - -t == =, - -0-65, 6<-, (7.8) 
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Figure 4: The coefficient K l=3 (Q) versus tan(9) in Ecliptic coordinates. Black line corresponds 
to H — 15°, the red line indicate the model with H = 30°, and the blue line corresponds to the 
h = 70°- KBOz model. 
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Figure 5: The coefficient K l (8) versus I. Black line corresponds to H — 15° and the red line 
indicate the model with H = 30°. 

7.2 Normalization on I = 5 component 

As alternative to the normalization of the direction of dipole modulation discussed above, 
we will in this section use a completely different approach, based on the minimization of 
power for the I = 5 multipole in the intrinsic CMB. The I = 5 mode of the power spectrum 
is one of the major sources of the parity asymmetry, clearly seen in Fig. || (see next 
section). As in the previous section, we will exploit the fact that only the £50 and e§ \ 
components of the modulation are non-zero. Taking under consideration, that the phase 
of the 5, 1 component of the ILC 7 is cp^i = —2.23845, we have adopted $ = —995,1 for 
the azimuthal angle $>, maximizing the correlations between the 5, 1 component of the ILC 
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Figure 6: Left column. Reconstructed ILC for the model with H = 30° (in Ecliptic coordinates, 
top panel). Second from the top is the same signal as the first, but in the Galactic coordinates. The 
second from the bottom panel, is the map of ILC modulated by the KBO (in Ecliptic coordinates). 
The bottom panel is the same map as above, just in Galactic coordinates. All the maps corresponds 
to normalization = tt, $ = —1.227 rad. Right column. The same as left, but for normalization 
on the I = 5 harmonic with 8 = it/2, $ = 2.2384 rad. 



and £5^. Then, from Eq.(^J) we get: 
K-(l = 5)(Q) = 



o^o/UcosG + 2|a5 i|^sinO 



11C (I = 5) [n 2 cos 2 9 + 2u 2 sin 2 9] 



(7.10) 
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Figure 7 : Left column. The map from the top left of Fig. |g, but for / = 2, 3, 4, 5 (from the top to 
the bottom). Right column. The map from the top right of Fig. but for Z = 2,3,4,5 (from the 
top to the bottom). 



where £5^ = /icosG and | ^5, 1 1 = z^sinG (see Eq.(7.4) and C(l = 5) is the power of ILC for 
the I = 5 mode. From Eq.( [7.10 ) one finds that the maximum value of K~(l = 5)(G) — 0.28 
can be achieved, if ~ ir/2. 
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Table 3: as im coefficients of the ILC 7 in the Ecliptic coordinates. 



8. De-correlation of ILC and KBO 



In this section we would like to discuss the possible changes in the morphology of the 
ILC map after removal of the modulations, associated with the KBO emission. We will 
concentrate on one particular model of the KBO emission, H = 30°, using two models of 
normalization for O and <£, discussed in the previous section. Hereafter we will call them 
model-13 and model-15. 

Firstly, we would like to point out that our de-correlation technique is generally unstable, 
since we cannot recover all the chance-correlations of intrinsic CMB and the KBO, which 
are sufficient for the low multipole range of the power spectrum (see |5(| for details). 
This is why de-correlation of the ILC and KBO could only demonstrate the tendency of 
the changes in the morphology of the ILC, pointing at some general properties of the 
reconstructed (de-correlated) signal. 

Secondly, it would be very naive to expect that the simple model of the KBO emission, 
which is based on a very general symmetry of KBO distribution in space, can significantly 
improve the ILC map for wide range of multipoles. However, for some particular multipoles 
even our simple model discussed in Section 2, could be very informative with respect to 
the possible direction of change in the morphology of the ILC map. 

Let us discuss this issue in greater detail. Our first step is to assume that the intrinsic CMB, 
associated with de-correlated ILC map, has zero correlations with the KBO emissivity map, 
which allows us to set k(1) = in Eq.(5.5). Then, taking Eq.fl4.2j) into account, for the 
even and odd components independently, we can get the following equation for c(Z) and 
p(l): 



*(0 = 



Er 



a l,rnXl m + a LmXW 



2 J2m \Xl,r, 



-K+(l) 



Lm\ 



p(iy- 



' a Lm £ l,m 



-K 



(0( 



2EmN,m| 2 VEr 

Then, taking the particular models for Xl,m an d £j m> discussed in the previous section, we 
can estimate the reconstructed signal from ILC 7, shown in Fig. |6| and Fig. ^ for model-13 
and model-15. In Fig. [8| we have plotted the power spectrum of the reconstructed signal 
with respect to the ILC 7 power. 
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Figure 8: Power spectrum for the reconstructed CMB in the model with H — 30° and normaliza- 
tion of the octupole (the red line) , and for normalization of I — 5 (the blue stars) and for ILC 7 
(the black line). 



9. Alignment and parity tests for de- correlated CMB 



We now perform a test of the alignment between the quadrupole and octupole for the de- 
correlated CMB signals, using the publicly available code from {http : / / www.phys.cwru.edu/projects/mpvector 
for the multipole vector approach from [51|. Our goal is to test, what the orientation is 
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Table 4: ai. m and a^. m coefficients of de-corrclatcd ILC 7 in Galactic coordinates for the models 
with normalization 13 and 15. 

between the various quadrupole from Table |I] and the octupole, and in particular the dif- 
ference between the alignment for pure ILC 7 and for our KBO-foreground cleaned maps. 
In [ 51 ] , the authors introduce various approaches for comparing the vectors associated with 
two different multipoles. We use the statistic of 'oriented area', following the definition from 
1 51], as summarized below. 



(9.1) 



where u( J1 >*) and v^J) come from the 11 multipole and v^ l2 ' k ^ and £>(' 2 ' m ) come from the 12 
multipole. That is, we cross the i'th. and j'th vector from multipole 11, and the fc'th and 
m'th vector from multipole 12, before we finally take the dot product between these two 
surfaces. Values close to 1 symbolizes a high level of alignment between the two planes. 
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Generally for a given 11 and 12 and % 7^ j and k 7^ m, there are M = /1(Z1 — l)/2(/2 — l)/4 
different products, meaning that for the comparison of quadrupole and octupole, we have 
three oriented areas (M1,M2 and M3). 

In Table ||, we have summarized the result of the quadrupole-octupole alignment test for the 
de-correlated CMB signal, including the standard result for the ILC 7 map for comparison. 
As is evident from Table ||, we have the usual strong alignment between quadrupole and 
octupole for the ILC. On the other hand, for the KBO-foreground cleaned signal, we see 
significantly reduced alignments among the three oriented areas. In summary, when we 
remove the KBO contribution from the ILC signal, to reduce the parity asymmetry, we see 
a weaker (statistically negligible) correlation between the quadrupole and octupole, than is 
the case for the ILC 7 map. Finally, for these two models (13 and 15) we have estimated 
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0.4433 
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Table 5: Result of the oriented area test between the quadrupole- and octupole, for the ILC case 
(a;,m) and the KBO cleaned signal (c; !m ), with 13 and 15 normalization from Section 7. 
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Figure 9: P-value for reconstructed CMB for the model with H = 30° and normalization on the 
octupole ^3(the green line), for normalization on / = 5 (the red line) and for ILC 7 (the blue line). 

the p-value to obtain a parity asymmetric power spectrum from pure random Gaussian 
and statistically isotropic CMB maps, shown in Fig. |9[ As one can see from this figure, 
de-correlation of the ILC 7 by the dipole modulated KBO can significantly increases the 
p-value by more than one order of magnitude for Z max ~ 20 — 25. At the end of this section 
we would like to point out that the amplitude of the dipole modulation for both these 
normalizations can be estimated from Fig. ^, with A ~ 0.3 — 0.4, which is comparable with 
the results by 0] at the < 2cr-level. 
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10. Conclusion 



In this paper, we have investigated a possible solution to the anomalous parity asymme- 
try in the Cosmic Microwave Background and the quadrupole-octupole alignment, related 
to the solar system foreground. The Kuiper Belt objects in the solar system plane can 
contribute to the microwave sky, and thus create a power contrast between even and odd 
multipoles. We have built a model of the KBO emissivity, based on the symmetry of the 
KBO in ecliptic coordinates. An essential part of the model is the incorporation of the 
dipole modulation of the ILC 7 map and KBO foreground, which transform the symmetric 
KBO foreground into an asymmetric part. An important point of our analysis is that by 
maximizing the cross-correlations between that asymmetric component and the ILC 7 odd 
harmonics for I = 3, and separately for Z = 5, we can fix the direction of dipole modulation, 
and significantly change the balance between even and odd multipoles in the ILC map. 
To illustrate the possible changes to the morphology of the low multipole domain of the 
CMB power spectrum, we have applied the strongest de-correlations between the ILC 7 
map and dipole modulated KBO foreground, in order to illuminate the clear improvement 
of the parity asymmetry and the quadrupole-octupole alignment. 

We have developed a method for cleaning the ILC 7 a^ m -values from the contribution of 
the dipole modulated KBO, thus retaining only the intrinsic CMB signal in the Ecliptic 
plane. It is possible to apply different weights to this filtering, allowing for various levels 
of parity symmetry for the low multipoles, depending on the normalization of the dipole 
modulation direction in the sky, and the corresponding amplitude of modulation. We have 
shown, that at the level where the parity is effectively restored, the quadrupole-octupole 
alignment would be significantly reduced down to the level of a chance correlation. 
In conclusion, the removal of the KBO contribution, in the framework of the dipole modula- 
tion model, clearly requires a statistical anisotropy of the CMB of non-cosmological origin, 
at least at the level of the dipole component. We believe, that the coming PLANCK data, 
with different systematics compared to the WMAP experiment, can put a light on the 
problem of low multipole anomalies in the CMB. 
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